NOLTR  64-163 


NOTICE 


Requests  for  additional  copies  by  Agencies  of  the  Department  of  Defense, 
their  contractors,  an;d  other  Government  agencies  should  be  directed  to: 

Defense  Documentation  Center  (DDC) 

Cameron  Station 
Alexandria,  Virginia 

Department  of  Defense  contractors1  who  have  established  DDC  services  or 
have  their  •need-to-know'  certified  by  the  coguizant  military  agency  of 
their  project  or  contract  should  ^lso  request  copies  from  DDC. 

All  other  persons  and  organizations  should  apply  to: 

U.  S.  Department  of  Commerce 
Office  of  Technical  Services 
Washington  25,  D.  C. 


NOLTR  64-163 


* 


MEASUREMENTS  AND  CORRELATION  OF  HEAT  TRANSFER 
IN  A  SOLID  PROPELLANT  ROCKET  NOZZLE 
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ABSTRACT:  The  two-dimensional  transient  heat  flow  in  a  conical 
De  Laval  nozzle  subject  to  a  solid  propellant  exhaust  flow  was 
investigated.  Local  temperature  and  heat  transfer  rates  at  the 
internal  surface  were  determined  from  a  combined  experimental 
and  analytical  method  developed  at  the  U.  S.  Naval  Ordnance 
Laboratory.  It  was  found  that  an  analysis  assuming  one¬ 
dimensional  heat  flow  gave  essentially  the  same  results  as  the 
two-dimensional  analysis. 

Local  heat  transfer  coefficients  were  compared  with  the  predic¬ 
tions  from  a  detailed  solution  of  the  turbulent  boundary  layer 
flow  and  with  the  widely  used  empirical  relation  of  Bartz.  The 
present  results  were  in  better  agreement  with  the  boundary  layer 
solution,  particularly  in  the  prediction  of  the  peak  value  at 
the  nozzle  throat. 

The  present  rocket  nozzle  data  expressed  in  Nusselt  number  form 
showed  very  good  agreement  with  the  steady  state  heat  transfer 
data  observed  with  a  water-cooled  copper  nozzle.  The  copper 
nozzle  data  was  measured  at  the  U.  S.  Naval  Ordnance  Laboratory 
using  heated  compressed  air. 

On  the  basis  of  the  present  data  plus  the  copper  nozzle  data  a 
simple  adjustment  is  derived  to  include  the  effects  of  compres¬ 
sibility  and  geometry  in  the  conventional  Nusselt-Reynolds 
number  correlation.  This  adjustment  decreases  the  scatter  of 
experimental  data  by  43  percent. 
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SYMBOLS 

A  cross  section  flow  area  of  nozzle 

Cp  specific  heat  of  gas  at  constant  pressure 

Cf  local  skin  friction  coefficient,  2^/p^u^3 

Cf.  incompressible  skin  friction  coefficient  for  zero 

1  heat  transfer  based  on  free  stream  conditions 

d  diameter 

h  heat  transfer  coefficient,  q/(Taw-Ts) 

H  enthalpy 

H'  reference  enthalpy 

k  thermal  conductivity 

M  Mach  number 

Nu^  Nusselt  number,  hd/k^ 

Nu  effective  Nusselt  number,  see  equation  (11) 

Nu^  incompressible  Nusselt  number,  see  equation  (9) 

Num  modified  Nusselt  number,  see  equation  (10) 

q  time  rate  of  heat  transfer  per  unit  area 

Pr  Prandtl  number,  ca  4  /k 

r  <x>  oo  00 

r  coordinate  in  radial  direction 

Reo  Reynolds  number,  p^u^d/p,* 

Stw  Stanton  number,  h/p^u^Cp^ 

Stj^  incompressible  Stanton  number,  i^c^Pr^"3/3 

T  temperature 

tj  radial  distance  from  nozzle  surface  to  thermocouple 

location,  see  figure  2  and  Table  I 

u  velocity 

x  axial  distance,  see  figure  2  and  Table  I 
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a  thermal  diffusivity 

y  ratio  of  specific  heats 

p  viscosity 

p  density  of  gas 

t  time 

tw  shear  stress 

<b  mass  flow 

Superscripts 

n  time  reference 

'  evaluated  at  Eckert's  reference  enthalpy 

*  nozzle  throat 

Subscripts 

aw  adiabatic  wall  condition 

i  space  reference  in  the  axial  direction 

j  space  reference  in  the  radial  direction 

r  radial  direction 

s  nozzle  surface  condition 

x  axial  direction 

o  stagnation  condition 

®  local  free-stream  condition 
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INTRODUCTION 

It  is  generally  accepted  that  the  exhaust  jet  from  a  solid 
propellant  contains  nongaseous  states  even  for  propellants 
without  solid  additives.  The  effect  of  these  nongaseous  states 
on  heat  transfer  to  motor  components  is  not  yet  well  defined 
(see,  for  example,  ref*  (1)).  Because  of  the  lack  of  a  good 
understanding  of  rocket  nozzle  heat  transfer  the  designer  often 
relies  on  empirical  relations  which  have  been  developed  for 
liquid  or  gaseous  flows  to  obtain  the  necessary  heat  transfer 
estimates.  It  follows  that  the  extension  of  these  empirical 
equations  to  determine  the  performance  of  rocket  motors  should 
be  verified  by  experimental  results.  It  can  be  seen  in  refer¬ 
ences  (2),  (3),  and  (4)  that  the  agreement  among  empirical 
relations  themselves  and  between  the  empirical  relations  and  the 
limited  experimental  data  available  can  vary  considerably 
depending  upon  the  choice  of  correlation  parameters.  These 
empirical  relations  can  be  improved  upon  or  at  least  the  limits 
of  application  can  be  determined  by  comparing  them  with  the 
results  from  carefully  controlled  experiments.  This  is  the 
objective  of  the  present  report. 

The  very  high  exhaust  temperatures  of  present-day  rocket 
motors,  greater  than  2000°K,  increase  the  difficulty  of 
measuring  heat  transfer  by  conventional  experimental  techniques. 
The  author  had  previously  developed  (ref.  (4))  a  combined 
experimental  and  analytical  method  for  obtaining  the  transient 
surface  temperature  and  heat  transfer  rate  at  the  throat  of  a 
solid  propellant  rocket  nozzle  where  heating  is  most  severe,. 

The  previous  method  assumed  one-dimensional  heat  flow  in  the 
radial  direction.  The  present  report  is  an  extension  of  the 
method  to  heat  flow  in  both  the  axial  and  radial  directions  to 
determine  the  surface  temperatures  and  heat  transfer  rates  over 
all  the  nozzle  surface.  In  an  attempt  to  correlate  the  data, 
the  free  stream  Nusseit  number  has  been  modified  for  compres¬ 
sibility  and  geometry. 


ANALYSIS 

A  combined  experimental-analytical  method  which  employs  a 
finite  difference  approximation  together  with  thermocouple  data 
was  developed  and  successfully  applied  in  reference  (4)  to 
describe  the  one-dimensional  heat  flow  at  the  throat  region  of 
a  solid  propellant  rocket  nozzle.  In  the  present  work  the 
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method  is  extended  to  describe  the  two-dimensional  heat  flow 
throughout  the  nozzle.  The  resulting  temperatures  and  heat 
transfer  rates  from  the  two-dimensional  solution  are  then 
compared  with  the  original  one-dimensional  solution  which 
assumes  heat  flow  only  in  the  radial  direction. 

Figure  1  illustrates  the  mesh  used  in  applying  the  f  ve 
difference  approximation.  The  region  R  is  defined  by  the 
location  of  thermocouples  on  its  boundary.  The  grid  lines  are 
orthogonal,  with  each  line  crossing  one  or  more  thermocouple 
locations.  Therefore,  the  grid  spacings  are  not  necessarily 
equal  but  are  determined  by  the  thermocouple  locations. 

The  thermocouples,  as  can  be  seen  in  figure  2  and  Table  1, 
have  been  mounted  no  closer  than  0.1  inch  from  the  inner  nozzle 
surface.  This  was  done  to  prevent  the  high  surface  temperatures 
from  destroying  the  thermocouples. 

The  heat  flow  within  region  R  is  described  by  the  two- 
dimensional  transient  equation  in  cylindrical  coordinates 
assuming  negligible  heat  flow  in  the  circumferential  direction: 


i  a 
r'3r 


„  n  dT 
cpP 


(1) 


Equation  (1)  is  solved  numerically  by  the  implicit  method  of 
finite  differences  (ref.  (5))  which  has  the  important  feature 
that  the  solution  is  stable  regardless  of  the  choice  of  time 
and  space  increments.  This  is  in  contrast  to  the  explicit 
method  which  is  widely  used  in  heat  flow  problems  and  which 
requires  a  stability  criterion  dependent  on  the  time  step  and 
mesh  spacing. 

Several  numerical  methods  are  available  for  the  solution 
of  two-dimensional  problems  in  the  implicit  form.  Some  of  these 
are  discussed  in  references  (6)  and  (7) .  The  method  selected 
is  an  alternating  direction  procedure  similar  to  that  presented 
in  reference  (8)  for  steady-state  problems  described  by 
elliptical  equations.  This  alternating  direction  procedure  may 
be  illustrated  as  follows.  (The  actual  grid  used  in  the  data 
reduction  is  rather  complicated,  but  the  main  ideas  may  be 
illustrated  with  the  simple  rectangular  mesh  shown  on  the  fol¬ 
lowing  page.) 
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Assuming  that  all  temperatures  at  time  interval  n  are  known,  and 
for  the  first  half  time-step  the  heat  conduction  in  the  r  direc¬ 
tion  is  neglected,  the  following  finite  difference  equation  in 
explicit  form  can  be  written  for  each  of  the  grid  points  indicat¬ 
ed  by  the  crosses:  (This  neglect  of  the  r  heat  conduction  is 
equivalent  to  dropping  the  r  derivatives  in  equation  (1). 


This  results  in  12  equations  in  terms  of  20  temperatures.  The 
eight  temperatures  indicated  by  circles  are  known  from  the 
thermocouple  readings.  Hence,  the  12  unknown  temperatures  at 
the  time  interval  (n+£)  and  at  the  position  indicated  by  the 
crosses  may  be  found  from  the  12  equations. 

For  the  next  half  time-step  at  time  interval  (n+1) ,  the 
procedure  is  repeated  but  a  somewhat  different  set  of  mesh  points 
is  used  in  formulating  the  equations.  The  mesh  points  are 
indicated  by  crosses,  and  the  heat  conduction  in  the  x  direction 
is  neglected. 
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The  finite  difference  equation 
for  each  of  the  grid  points  is 


1 

1  1 

rj 

k. 

n+£ 

n+i 

k  r 

(j-i) 

lT(j-o  -  Tj 

rj  "  r(j-0 


for  the  second  half  time  interval 
then: 


(3) 


c+£ 


MC 


n+i  n+£ 

V  _  T 

LhJ—Jul 

At 


Similar  to  before,  the  temperature  at  each  of  the  ten  stations 
indicated  by  a  cross  at  time  interval  (n+1)  may  be  computed  from 
equation  (3),  At  time  interval  (n+l|),  the  procedure  is  similar 
to  time  interval  (n4-|) .  The  alternating  procedure  is  repeated 
at  all  succeeding  times. 

This  alternating  procedure  which  is  equivalent  to  computing 
the  one-dimensional  heat  flow  in  alternate  directions  simplifies 
the  machine  calculation  on  the  IBM  7090.  The  numerical  problem 
is  reducea  to  the  solution  of  a  tri-diagonal  matrix  which  can 
be  readily  solved  by  Gaussian  elimination. 

The  temperature  at  the  inner  nozzle  surface  is  extrapolated 
from  the  known  results  in  region  R  by  an  equation  similar  to 
equation  (3) .  The  use  of  equation  (3)  for  extrapolation  will 
introduce  errors  which  can  be  kept  to  a  minimum  by  maintaining 
the  extrapolation  distance  small.  The  machine  program  selects 
a  grid  at  each  extrapolation  location.  The  surface  heat  transfer 
rates  are  determined  from  the  temperature  gradient  at  the 
surface . 


EXPERIMENTAL  PROCEDURE  AND  INSTRUMENTATION 

Rocket  nozzle  firings  were  conducted  at  the  Rocket  Tunnel 
Facility  of  the  Applied  Physics  Laboratory,  Johns  Hopkins 
University.  The  propellant  used  was  of  the  standard  double  base, 
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end  burning,  10-second  grain  supplied  by  the  Allegany  Ballistics 
Laboratory.  This  grain  produces  a  nominal  operating  supply 
pressure  and  temperature  of  1150  psia  and  2500°K,  respectively, 
in  the  combustion  chamber.  A  detailed  discussion  of  the  trans¬ 
port  and  thermodynamic  properties  of  this  propellant  is  given 
in  references  (9)  and  (10) . 

A  diagram  of  the  heat  transfer  model  consisting  of  a  solid, 
heavy-wall  molybdenum  heat  sink  nozzle  is  shown  in  figure  2. 

The  internal  nozzle  geometry  is  the  standard  conical  nozzle 
configuration  used  at  the  APL  facility.  The  pertinent  dimensions 
are  given  in  Table  1  and  shown  in  figure  2. 

The  nozzle  was  instrumented  with  a  total  of  forty-five 
thermocouples  placed  in  nine  axial  locations  as  shown.  The 
peripheral  thermocouples  were  used  to  define  the  boundary  for 
the  region  R  described  previously  and  also  the  mesh  spacing  for 
the  two-dimensional  numerical  solution.  The  internal  thermo¬ 
couples  were  used  for  checking  the  results  of  the  one  and  two- 
dimensional  numerical  solution. 

All  thermocouples  were  made  from  platinum  and  platinum- 
rhodium  30  gage  wires.  The  junctions  were  spot-welded  onto 
one-degree  tapered  molybdenum  plugs,  with  one  plug  for  each 
axial  location.  Each  plug  was  inserted  into  a  hand-fitted 
mating  hole  which  bottomed  at  approximately  .01  inch  from  the 
internal  surface.  All  the  thermocouples  were  located  in  one 
plane  that  passed  through  the  nozzle  axis.  The  installation  of 
the  wires  into  the  tapered  plugs  was  similar  to  that  reported 
in  reference  (4) .  The  emf  developed  by  the  thermocouples  were 
recorded  on  two  50-Channel  Midwestern  Direct  Recording  Oscillo¬ 
graphs. 

The  supply  pressure  was  measured  with  a  pressure  transducer 
connected  to  a  static  orifice  located  at  the  downstream  end  of 
the  combustion  chamber.  The  supply  temperature  was  measured 
with  an  unshielded  .020-inch  diameter  tungsten-iridium  thermo¬ 
couple  located  in  the  flow  directly  ahead  of  the  nozzle  inlet. 

The  nozzle  was  attached  to  the  combustion  chamber  by  a 
slotted-thread  connector  wh.ch  needed  only  a  quarter  turn  of  the 
nozzle  to  seal  the  passage.  This  quick  connector  simplified 
the  installation  of  the  nozzle  and  prevented  damage  to  the 
delicate  thermocouple  wires. 


DISCUSSION  OF  DATA 

The  characteristic  pressure  and  temperature  histories  in 
the  combustion  chamber  are  shown  in  figures  3  and  4,  respectively. 
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The  various  "dips"  in  the  temperature  time  curve  of  figure  4, 
as  concluded  from  previous  tests  (see  ref.  (4)),  are  caused  by 
the  temporary  insulation  effect  of  nongaseous  products  deposited 
on  the  bare  thermocouple.  For  the  heat  transfer  analysis,  the 
dashed  curve  through  the  maximum  values  represents  the  supply 
temperature  of  the  flow. 

Continuous  temperature-versus-time  recordings  of  each  of 
the  forty-five  thermocouples  were  taken  during  the  course  of 
the  9.5-second  run.  For  the  numerical  solution  the  data  at 
time  increments  of  0.1  and  0.25  seconds  were  used.  The  dif¬ 
ference  in  temperature  distribution  as  calculated  with  these 
two  time  increments  is  negligible. 

The  correlation  between  the  experimental  results  and  the 
one  and  two-dimensional  numerical  methods  were  made  at  three 
selected  time  stations.  These  were  at  three  seconds;  when 
free-stream  flow  conditions  should  have  stabilized;  at  an 
intermediate  station  of  six  seconds;  and  just  before  burnout  at 
9.5  seconds.  The  temperature  distributions  at  each  of  the  nine 
axial  locations  at  these  three  time  stations  are  shown  in 
figures  5  through  13.  The  computed  two-dimensional  solution 
showed  good  agreement  with  measured  data  at  the  9,5-second  time 
station.  At  the  early  times  of  three  and  six  seconds,  the 
computed  two-dimensional  solution  appeared  to  be  lower  than 
experimental  data;  the  lag  was  more  severe  at  three  seconds, 
particularly  at  the  thicker  wall  stations.  Since  there  is  no 
reason  to  doubt  the  experimental  data,  the  explanation  for  this 
discrepancy  is  that  the  alternating  procedure  selected  for 
solving  equation  (2)  introduces  a  time  lag  in  the  computed  tem¬ 
peratures  at  the  early  highly  transient  conditions. 

Figure  14  is  a  graph  of  the  computed  surface  temperature 
at  times  3  and  9.5  seconds.  In  addition  to  the  anticipated 
high  temperature  at  the  nozzle  throat  area,  a  second  hot  spot 
occurred  at  slightly  less  than  one  inch  downstream  from  the 
nozzle  throat  (station  G) .  This  second  temperature  maximum  can 
also  be  seen  in  the  map  of  the  isotherms  at  6  and  9.5  seconds 
shown  in  figures  15  and  16,  respectively.  The  exact  cause  of 
this  second  maximum  is  not  known;  however,  it  may  be  associated 
with  the  flow  separation  and  reattachment  observed  with  the 
same  nozzle  geometry  in  other  tests. 

It  is  noted  in  figures  15  and  16  that  the  calculated  iso¬ 
therms  are  wavy  near  the  inner  surface.  As  one  goes  away  from 
this  surface,  the  waviness  is  retained  in  the  one-dimensional 
curves  but  is  diminished  in  the  two-dimensional  ones.  One 
expects  the  waviness  to  diminish  as  one  gets  further  from  the 
hot  surface — hence.,  it  appears  that  the  two-dimensional  calcula¬ 
tion  is  a  more  accurate  description  of  the  heat  conduction  process 
than  the  one-dimensional. 
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Local  nozzle  surface  heat  transfer  coefficients  are  shown 
in  figure  17.  The  rocket  nozzle  data  show  fair  agreement  with 
the  turbulent  boundary  layer  solutions  of  Persh  and  Lee, 

(ref.  (11))  which  assumed  ideal  gas  flow  and  constant  wall  tem¬ 
peratures  of  555°K  and  2222°K.  Also  plotted  is  the  widely  used 
solution  of  Bartz  (ref.  (12))  which  gives  a  higher  heat  transfer 
coefficient  at  the  throat  and  downstream.  Unlike  the  ideal  gas 
flow  calculation  of  reference  (11)  where  maximum  heat  transfer 
is  predicted  to  occur  slightly  upstream  of  the  nozzle  throat, 
the  rocket  nozzle  results  show  the  peak  to  occur  slightly  down¬ 
stream.  The  upstream  heat  transfer  coefficients  are  lower  than 
the  boundary  layer  prediction  possibly  due  to  the  observed  large 
combustion  deposits  in  this  region.  The  high  data  point  at 
approximately  0.9  inch  is  due  to  the  probable  reattachment  of 
separated  flow  discussed  earlier. 

For  engineering  purposes,  the  maximum  temperature  and  the 
maximum  stress  are  of  prime  importance.  Both  of  these  occur  at 
the  nozzle  surface.  The  computed  results  show  that  the  difference 
between  the  one-dimensional  solution  and  the  two-dimensional 
solution  in  computing  the  surface  temperature  and  heat  transfer 
are  within  experimental  accuracy  for  the  test  configuration  used. 


CORRELATION  OF  DATA 

The  conventional  correlation  of  nozzle  heat  transfer  data 
in  terms  of  Nusselt  number  versus  Reynolds  number  is  shown  in 
figure  18.  The  steady  state  heat  transfer  data  of  reference  (13), 
which  covered  an  extensive  range  of  Reynolds  number,  is  also 
plotted  and  appeared  to  be  in  very  good  agreement  with  the 
present  rocket  nozzle  data.  The  scatter  of  the  data  may  be 
represented  by  the  distance  between  the  two  dashed  lines  drawn. 
The  upper  line  is  drawn  through  the  higher  data  points  and  the 
lower  line  is  drawn  through  the  lowest  data  point  parallel  to 
the  upper  line.  The  two  dashed  lines  deviate  from  their  arith¬ 
metic  mean  by  +  49  percent. 

It  has  been  noted  in  reference  (13)  for  air  flow  and  also 
reported  in  reference  (2)  for  a  simulated  liquid  propellant  flow 
that  the  heat  transfer  rates  are  higher  in  the  subsonic  nozzle 
inlet  region  than  at  the  supersonic  expansion  region  for  the 
same  Reynolds  number.  This  is  also  true  for  the  present  rocket 
nozzle  data  and  is  illustrated  in  figure  18.  The  subsonic  data 
are  shown  by  the  hollow  symbols  for  both  the  air  flow  data  of 
reference  (13)  and  for  the  present  rocket  nozzle  flow  data. 

These  appear  to  be  in  good  agreement  with  the  subsonic  turbulent 
pipe  flow  relation  of  Dittus  and  Boelter  (ref.  (14))  which  is 
represented  by  tho  line  of  the  following  expression: 
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Nu  »  0.0265  Re0-8  Pr0-3  (4) 

From  the  observations  described  in  the  preceding  paragraph, 
it  appears  that  a  compressibility  correction  may  reduce  the 
scatter  of  the  data.  This  compressibility  correction  is  intended 
as  an  attempt  to  remove  the  Mach  number  and  wall  temperature 
effects  but  not  the  experimental  scatter  or  the  difference 
between  the  one-  and  two-dimensional  solutions  or  between  solu¬ 
tions  for  different  times. 

Two  relatively  simple  compressibility  correction  methods 
which  have  been  successfully  applied  to  flat  plate  flow  were 
tried  here  for  nozzle  flow.  The  first  is  by  evaluating  the  data 
at  a  reference  enthalpy  (see  ref.  (15))  which  is  defined  as: 


-  1  (5) 


The  results  of  applying  the  reference  enthalpy  method, 
shown  in  figure  19,  showed  no  improvement  in  the  scatter  of  the 
rocket  nozzle  data.  However,  the  cold-flow  data  of  reference 
(13)  was  reduced  from  the  original  scatter  of  +  44.6  percent 
when  correlated  in  terms  of  free-stream  properties  to  a  scatter 
of  +  35.7  percent  when  correlated  by  the  reference  enthalpy 
metEod.  A  possible  reason  for  the  lack  of  favorable  results 
for  the  rocket  nozzle  data  is  that  the  constants  in  equation 
(5)  were  empirically  determined  for  a  particular  set  of  condi¬ 
tions.  These  original  constants,  which  were  adjusted  to  fit 
low  temperature  flat  plate  data  as  demonstrated  in  reference  (16) , 
are  not  necessarily  correct  for  high  temperature  nozzle  flows  in 
which  a  pressure  gradient  exists.  However,  as  shown,  the  method 
does  reduce  somewhat  the  scatter  of  cold  flow  data  with  pressure 
gradient. 


~  =  1  +  0.032  Mi  +0.58 
Hoo 


The  second  method  is  by  the  empirical  Winkler-Cha  skin 
friction  formula  (see  ref.  (16))  where  the  local  skin  friction 
coefficient  can  be  approximated  by: 


(6) 


where  Cf ^  is  the  incompressible  skin  friction  coefficient,  which 
is  a  known  function  of  the  momentum  thickness  Reynolds  number. 
The  right  hand  side  of  equation  (6)  may  be  looked  upon  as  a 
function  of  Mach  number  and  wall  temperature  which  when  multi¬ 
plied  by  Cf  removes  the  Mach  number  and  wall  temperature 
dependencies.  If  one  further  assumes  that  Colburn's  version  of 
Reynolds  analogy  holds,  that  is 
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cf 


00 


2 


3/3 

Stjr. 


(7) 


then  it  may  be  shown  that 


Nui 


IZa]* 

Taw 

UJ 

i  TS 

where 

Nui  5  sti  ReJPr* 


(8) 


(9) 


In  other  words,  the  same  factors  which  remove  these  variations 
from  the  skin  friction  coefficient  for  a  flat  plate  should  also 
remove  them  from  the  Nusselt  number. 


The  results  of  using  the  Winkler-Cha  method  are  shown  in 
figure  20.  Similar  to  the  result  from  the  reference  enthalpy 
method,  the  rocket  nozzle  data  showed  negligible  change  while 
the  low  temperature  data  was  reduced  from  +  44.6  percent  to 
f  34.1  percent.  It  should  be  noted  that  both  the  reference 
enthalpy  and  the  Winkler-Cha  methods  give  nearly  the  same 
improvement  to  the  low  temperature  data  and  as  shown  in 
figures  19  and  20  put  the  low  temperature  data  in  line  with  the 
higher  data  points  of  the  rocket  nozzle  results. 

Another  method  of  correlation  was  carried  out  as  follows. i 
In  analogy  with  figure  22  of  reference  (17),  the  term  (Taw/Ts) 4 
was  omitted  from  equation  (6).  We  define 

Nu™  "(t*)4  Nu”  (10) 

Fig.  21  shows  the  results.  The  overall  scatter  of  both  the  low 
temperature  and  rocket  nozzle  data  in  figure  21  (excluding  the 
data  at  the  second  maximum  temperature  region  which  appeared 
outside  the  band  drawn)  was  reduced  from  the  original  +  49 
percent  to  36  percent.  Comparison  of  figures  18  through  21  shows 
that  of  the  four  methods  tried,  this  modified  version  of  the 
Winkler-Cha  formula,  equation  (10),  does  the  best  job  in  reducing 
the  scatter  in  the  heat  transfer  data. 

A  closer  look  at  the  nozzle  throat  data  of  figure  21,  which 
is  replotted  in  figure  22,  shows  that  a  straight  line  can  be 
drawn  through  these  data  points.  Further  analysis  of  the  sub¬ 
sonic  data  showed  a  systematically  increasing  deviation  from 
this  straight  line  with  increasing  distance  (decreasing  Reynolds 
number)  from  the  nozzle  throat.  This  is  illustrated  by  the 
crosses  on  figure  22  which  represent  the  data  of  reference  (13) 
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for  one  test  condition.  The  dashed  line  indicates  the  trend  of 
the  subsonic  data.  Since  the  free  stream  temperature  changed 
very  little  in  the  subsonic  region,  the  modified  Winkler-Cha 
compressibility  correction  is  negligible  at  the  nozzle  inlet. 

An  inspection  of  the  parameters  used  in  the  correlation  suggest 
that  the  nozzle  geometry  may  not  be  fully  accounted  for  in  the 
applied  correction.  Therefore,  an  added  geometry  correction 
factor  in  the  form  of  (d/d*)^  was  tried  and  appeared  to  work 
satisfactorily.  The  exponent  £  was  determined  by  calculating 
d/d*  and  adjusting  £  by  trial  until  the  data  points  fell  nearest 
to  the  line  drawn  through  the  nozzle  throat  data  points.  This 
resulted  in  a  value  of  t  equal  to  0.6.  Combining  the  above 
geometry  factor  together  with  the  modified  compressibility 
factor  of  Winkler-Cha  results  in  the  following  adjustment  in  the 
Nusselt  number  correlation: 

0.&  C.6 


Nu. 


Nu 


For  isentropic  one-dimensional  flows  where  the  diameter  is 
related  to  the  Mach  number,  equation  (11)  can  be  written  in  a 
different  form: 


Nu, 


Y+l 

2 


Y-l 

2 


u.2 

In- 


7Y—  13 

3TT  (y-1} 


(12) 


A  constant  effective  value  of  y  “  1.25  may  be  assumed.  This 
value  is  computed  at  the  mean  temperature  between  stagnation 
and  nozzle  throat  surface  temperatures  (see  refs.  (9)  and  (10)). 

The  graph  of  the  effective  Nusselt  number  versus  Reynolds 
number  is  shown  in  figure  23.  The  deviation  of  the  data  from 
the  mean  value  is  reduced  from  the  original  +  49  percent  to 
+  28  percent.  A  line  drawn  through  the  data  would  result  in 
The  following  expression: 

Nue  =  0.001  Re^  (13) 

The  foregoing  adjustment  of  conventional  correlation  para¬ 
meters  is  an  attempt  to  reduce  the  wide  scatter  of  rocket  nozzle 
heat  transfer  data.  The  extent  of  applicability  of  this  cor¬ 
relation  will  need  further  experimental  support.  Equations  (11) , 
(12),  and  (13)  were  derived  from  data  which  extended  over  the 
following  ranges: 


4.0  x  105  «:  Rem  <  4.1  x  106 

CD 

-1.0  *  (T0/TJ  ^  2.x 
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1.4  *  (Taw/Ts)  «.  10.3 
0.35  s  (d*/d)  *  1.0 

From  equations  (11)  and  (13)  and  the  definition  of  the 
Stanton  number,  the  following  equation  can  be  derived  to  compute 
the  local  heat  transfer  coefficient: 


In  many  instances  the  specific  heat  and  Prandtl  number  change 
verv  little  throughout  the  flow  temperature  ranu'e.  Equation 
(1  •  can  then  be  reduced  further  to  a  constant  plus  three  terms 
which  involve  only  the  mass  flow,  local  diameter,  and  the  local 
Mach  number.  The  local  heat  transfer  coefficient  computed  by 
equation  (14)  is  compared  with  the  present  rocket  nozzle  data 
in  figure  17. 


SUMMARY  AND  CONCLUSIONS 

Detailed  investigations  of  the  two-dimensional  transient 
heat  flow  have  been  made  for  a  heat  sink  nozzle  subject  to  a 
solid  propellant  flow.  A  finite  difference  implicit  numerical 
method  was  employed  to  describe  the  two-dimensional  heat  con¬ 
duction  and  to  obtain  the  internal  surface  heat  transfer  rates. 
The  surface  temperature  and  heat  transfer  were  compared  and 
were  found  to  agree  with  the  results  from  the  one-dimensional 
heat  flow  analysis  to  within  the  measuring  accuracy. 

Peak  nozzle  throat  heat  transfer  coefficients  agree  with 
the  predictions  from  a  detailed  turbulent  boundary  layer 
solution.  Propellant  deposits  and  flow  separation  and  reattach¬ 
ment  affected  the  heat  transfer  upstream  and  downstream  of  the 
nozzle  throat  respectively.  Nusselt  number  correlation  showed 
very  good  agreement  with  steady  state  heat  transfer  data 
observed  with  a  water-cooled  copper  nozzle. 

On  the  basis  of  the  present  data  plus  the  copper  nozzle 
data  a  simple  adjustment  is  derived  to  include  the  effects  of 
compressibility  and  geometry  in  the  conventional  Nusselt- 
Reynolds  number  correlation.  This  adjustment  decreases  the 
scatter  of  experimental  data  by  43  percent.  The  data  can  be 
approximated  by  the  expression 

Nup  =  0.001  Re  (13) 

C  CD 
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from  which  a  simple  equation  for  the  local  neat  transfer  coef¬ 
ficient  was  obtained. 


h 


=  0.001 

!T0] 

—  0.5 

-0.6 

1*11 

CPoo 

! 

iTwj 

i  d  I 

1a  /  1 

l*J 

(14) 


The  above  expression  is  applicable  over  the  following  ranges: 

4.0  x  10s  <  Re  <  4.1  x  106 

oo 

1.0  <:  (T0/T J  <;  2.1 
1.4  <  (Taw/Ts)  £  10.3 
0.35  s  (d*/d)  <:  1.0 
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TABLE  I 

Thermocouple  Locations 
(See  figure  2) 


Stations 

x,  inches 

tj,  inches 

A 

1.109 

.187 

.272 

.492 

972 

1.456 

B 

1.475 

.177 

.262 

.482 

1.162 

1.821 

C 

1.775 

.135 

.220 

.440 

1.120 

2.004 

D 

2.015 

.139 

.224 

.444 

1.124 

2.043 

E 

2.255 

.122 

.207 

.427 

1.107 

2.026 

F 

2.500 

.103 

.188 

.408 

1.088 

1.972 

G 

2.890 

.117 

.202 

.422 

1.102 

1.886 

H 

3.395 

.130 

.215 

.435 

1.115 

1.774 

I 

4.000 

.170 

.255 

.475 

.955 

1.639 

5. 154 
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FIG  2  SCHEMATIC  OF  INSTRUMENTED  ROCKET  NOZZLE 
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FIG  3  MEASURED  CHAMBER  PRESSURE  AS  A  FUNCTION  OF 
TIME  FROM  IGNITION 
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FIG  4  MEASURED  TOTAL  TEMPERATURE  AS  A  FUNCTION 
OF  TIME  FROM  IGNITION 
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FIG  5  TEMPERATURE  DISTRIBUTIONS  AT  STATION  A 
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FIG  6  TEMPERATURE  DISTRIBUTIONS  AT  STATION  B 


TEMPERATURE  (°K) 


NOLTR  64-163 


TEMPERATURE  (°K) 


NOLTR  64-163 


DISTANCE  FROM  NOZZLE  SURFACE  (INCHES) 

FIG  9  TEMPERATURE  DISTRIBUTIONS  AT  STATION  E 
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DISTANCE  FROM  NOZZLE  SURFACE  (INCHES) 
FIG  10  TEMPERATURE  DISTRIBUTIONS  AT  STATION  F 
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FIG  11  TEMPERATURE  DISTRIBUTIONS  AT  STATION  G 
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FIG  13  TEMPERATURE  DISTRIBUTIONS  AT  STATIONS  I 


FIG  14  NOZZLE  SURFACE  TEMPERATURE  DISTRIBUTIONS 
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FIG  16  NOZZLE  ISOTHERMS  9.5  SECONDS  AFTER  IGNITION 
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FIG  20  NUSSELT- REYNOLDS  DUMBER  CORRELATION  BASED  ON  FREE -STREAM 
PROPERTIES  AND  USING  THE  WINKLER- CHA  COMPRESSIBILITY  CORRECTION 
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Headquarters ,  Arnold  Engineering  Development  Center 
ARO,  Inc. 

Arnold  Air  Force  Station,  Tennessee 
Attn:  Technical  Library 
Attn:  AEOR 
Attn:  AEOIM 

Commanding  Officer,  Harry  Diamond  Laboratories 
Washington  25,  D.  C. 

Attn:  Library,  Room  211,  Bldg.  92 

Commanding  General 

U.  S.  Army  Missile  Command 

Redstone  Arsenal,  Alabama 

Attn:  AMSMI-RR  (Mr.  N.  Shapiro) 

Attn:  AMSMI-RB  (Redstone  Scientific  Information 

Center) 
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NASA 

George  C,  Marshall  Space  Flight  Center 
Huntsville,  Alabama 
Attn:  Dr.  E.  Geissler 
Attn:  Mr.  T.  Reed 
Attn:  Mr.  H.  Paul 
Attn:  Mr.  W.  Dahm 
Attn:  Mr.  H.  A.  Connell 
Attn:  Mr.  J.  Kingsbury 
Attn:  ARDAB-DA 

APL/JHU  (NOw  7386) 

8621  Georgia  Avenue 
Silver  Spring,  Maryland 

Attn:  Technical  Reports  Group  2 

Attn:  Mr.  Richard  Suess 

Attn:  Dr.  F.  Hill 

Attn:  Dr.  L.  L.  Cronvich 

Scientific  &  Technical  Information  Facility 
P.0.  Box  5700 
Bethesda,  Maryland 

Attn:  NASA  Representative  (S-AK/DL) 

Commander 

Air  Force  Flight  Test  Center 
Edwards  Air  Force  Ease 
Muroc,  California 
Attn :  FTOTL 

Air  Force  Office  of  Scientific  Research 
Hollomon  Air  Force  Base 
Alamogordo,  New  Mexico 
Attn:  SRLTL 

U.  S.  Army  Engineer  Research  &  Development 
Laboratories 
Fort  Belvoir,  Virginia 
Attn:  STINFO  Branch 
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University  of  Minnesota 
Minneapolis  14,  Minnesota 
Attn:  Dr.  E.  R.  G.  Eckert 
Attn:  Heat  Transfer  Laboratory 
Attn:  Technical  Library 

Rensselaer  Polytechnic  Institute 
Troy,  New  York 

Attn:  Dept,  of  Aeronautical  Engineering 

Dr.  James  P.  Hartnett 

Department  of  Mechanical  Engineering 

University  of  Delaware 

Newark,  Delaware 

Princeton  University 
James  Forrestal  Research  Center 
Gas  Dynamics  Laboratory 
Princeton,  New  Jersey 

Attn:  Prof.  S.  Bogdonoff 

Attn:  Dept,  of  Aeronautical  Engineering  Library 

Defense  Research  Laboratory 
The  University  of  Texas 
P.  0.  Box  8029 
Austin  12,  Texas 

Attn:  Assistant  Director 

Ohio  State  University 
Columbus  10,  Ohio 

Attn:  Security  Officer 

Attn :  Aerodynamics  Laboratory 

Attn:  Dr.  J.  Lee 

Attn:  Chairman,  Dept,  of  Aero.  Engineering 

California  Institute  of  Technology 
Pasadena,  California 

Attn:  Guggenheim  Aero.  Laboratory, 

Aeronautics  Library 
Attn:  Jet  Propulsion  Laboratory 
Attn:  Dr.  H.  Liepmann 
Attn:  Dr.  L.  Lees 
Attn:  Dr.  D.  Coles 
Attn:  Dr.  A.  Roshko 


Case  Institute  of  Technology 
Cleveland  6,  Ohio 
Attn:  G.  Kuerti 
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North  American  Aviation,  Inc. 

Aerophysics  Laboratory 
Downey,  California 

Attn:  Chief,  Aerophysics  Laboratory 
Attn:  Missile  Division  (Library) 

Department  of  Mechanical  Engineering 
Yale  University 
400  Temple  Street 
New  Haven,  Connecticut 
Attn:  Dr.  P.  P.  Wegener 

MIT  Lincoln  Laboratory 
Lexington,  Massachusetts 

RAND  Corporation 
1700  Main  Street 
Santa  Monica,  California 

Attn:  Library,  USAF  Project  RAND 

Attn:  Technical  Communications 

Mr.  J.  Lukasiewicz ,  Chief 
Gas  Dynamics  Facility 
ARO,  Incorporated 
Tullahoma,  Tennessee 

Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts 


Attn: 

Prof. 

J.  Kaye 

Attn: 

Prof. 

M,  Fins ton 

Attn : 

Mr.  J. 

Baron 

Attn : 

Prof. 

A.  H.  Shapiro 

Attn: 

Naval 

Supersonic  Laboratory 

Attn : 

Aero. 

Engineering  Library 

Attn: 

Prof. 

Ronald  F.  Probstein 

Attn : 

Prof . 

C,  C.  Lin 

Polytechnic  Institute  of  Brooklyn 
527  Atlantic  Avenue 
Freeport,  New  York 
Attn:  Dr.  M.  Bloom 

Attn:  Aerodynamics  Laboratory 


Brown  University 
Division  of  Engineering 
Providence,  Rhode  Island 
Attn:  Librarian 


2 


iv--4£iUU 


AERODYNAMICS  DEPARTMENT 
EXTERNAL  DISTRIBUTION  LIST  (A2) 

No.  of 
Copies 

Air  Ballistics  Laboratory 
Army  Ballistic  Missile  Agency 
Huntsville,  Alabama 

Applied  Mechanics  Reviews 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas 

BUWEPS  Representative 
Aerojet-General  Corporation 
6352  N.  Irwindale  Avenue 
Azusa,  California 

The  Boeing  Company 
Seattle,  Washington 

Attn:  J.  H.  Russell,  Aero-Space  Division 
Attn :  Research  Library 

United  Aircraft  Corporation 
400  Main  Street 
East  Hartford  8,  Connecticut 
Attn:  Chief  Librarian 

Attn:  Mr.  W.  Kuhrt,  Research  Department 

Attn:  Mr.  J.  G.  Lee  2 

Hughes  Aircraft  Company 
Florence  Avenue  at  Teale  Streets 
Culver  City,  California 
Attn:  Mr.  D.  J.  Johnson 

R&D  Technical  Library 

McDonnell  Aircraft  Corporation 

P,  0.  Box  516 

St.  Louis  3,  Missouri 

Lockheed  Missiles  and  Space  Company 
P.  0.  Box  504 
Sunnyvale,  California 
Attn:  Dr.  L.  H.  Wilson 
Attn:  Mr.  M.  Tucker 
Attn:  Dr.  R.  Smelt 

Martin  Company 
Baltimore,  Maryland 
Attn:  Library 
Attn:  Chief  Aerodynamic 1st 
Attn:  Dr,  W.  Morkovin,  Aerophysics  Division 
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CONVAIR 

A  Division  of  General  Dynamics  Corporation 
Fort  Worth,  Texas 
Attn:  Library 

Attn:  Theoretical  Aerodynamics  Group 
Purdue  University 

School  of  Aeronautical  &  Engineering  Sciences 
LaFayette,  Indiana 

Attn:  R.  L,  Taggart,  Library 

University  of  Maryland 
College  Park,  Maryland 

Attn:  Director  2 

Attn:  Dr.  J.  Burgers 

Attn:  Librarian,  Engr.  &  Physical  Sciences 
Attn:  Librarian,  Institute  for  Fluid  Dynamics 
and  Applied  Mathematics 
Attn:  Prof.  S.  I.  Pai 

University  of  Michigan 
Ann  Arbor,  Michigan 
Attn:  Dr.  a.  Kuethe 
Attn:  Dr.  At  Laporte 

Attn:  Department  of  Aeronautical  Engineering 

Stanford  University 
Palo  Alto,  California 

Attn:  Applied  Mathematics  &  Statistics  Lab. 

Attn:  Prof.  D.  Bershader,  Dept,  of  Aero.  Engr. 

Cornell  University 

Graduate  School  of  Aeronautical  Engineering 
Ithaca,  New  York 

Attn:  Prof.  W.  R.  Sears 

The  Johns  Hopkins  University 
Chrrles  and  34th  Streets 
Baltimore,  Maryland 

Attn:  Dr.  F.  H.  Clauser 

University  of  California 
Berkeley  4,  California 
Attn:  G.  Maslach 
Attn:  Dr.  S,  A.  Schaaf 
Attn:  Dr.  Holt 

Attn:  Institute  of  Engineering  Research 
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Cornell  Aeronautical  Laboratory,  Inc, 

4455  Genesee  Street 
Buffalo  21,  New  York 
Attn:  Librarian 
Attn:  Dr.  Franklin  Moore 
Attn:  Dr.  J.  G.  Hall 
Attn:  Mr.  A.  Hertzberg 

University  of  Minnesota 
Rosemount  Research  Laboratories 
Rosemount,  Minnesota 

Attn:  Technical  Library 

Director,  Air  University  Library 
Maxwell  Air  Force  Base,  Alabama 

Douglas  Aircraft  Company,  Inc. 

Santa  Monica  Division 
3000  Ocean  Park  Boulevard 
Santa  Monica,  California 

Attn:  Chief  Missiles  Engineer 
Attn:  Aerodynamics  Section 

CONVAIR 

A  Division  of  General  Dynamics  Corporation 
Daingerf ield,  Texas 

CONVAIR 

Scientific  Research  Laboratory 
5001  Kearney  Villa  Road 
San  Diego,  California 

Attn:  Asst,  to  the  Director  of  Scientific  Research 
Attn:  Dr.  B.  M.  Leadon 
Attn:  Library 

Republic  Aviation  Corporation 
Farmingdale,  New  York 

Attn:  Technical  Library 

General  Applied  Science  Laboratories,  Inc. 

Merrick  and  Stewart  Avenues 
Westbury,  L.  I.,  New  York 
Attn:  Mr.  Walter  Daskin 
Attn:  Mr.  R.  W.  Byrne 
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Arnold  Research  Organization,  Inc. 
Tullahoma,  Tennessee 

Attn:  Technical  Library 

Attn:  Chief,  Propulsion  Wind  Tunnel 

Attn:  Dr.  J.  L.  Potter 

General  Electric  Company 
Missile  Space  Division 
3198  Chestnut  Street 
Philadelphia,  Pennsylvania 

Attn:  Larry  Chasen,  Mgr.  Library 

Attn:  Mr.  R.  Kirby 

Attn:  Dr.  J.  Farber 

Attn:  Dr.  G.  Sutton 

Attn:  Dr.  J.  D.  Stewart 

Attn:  Dr.  S.  M.  Scala 

Attn:  Dr.  H.  Lew 

Attn:  Mr.  J.  Persh 


Eastman  Kodak  Company 
Navy  Ordnance  Division 
50  West  Main  Street 
Rochester  14,  New  York 
Attn:  W.  B.  Forman 

Library 

AVCO-Everett  Research  Laboratory 
2385  Revere  Beach  Parkway 
Everett  49,  Massachusetts 

Chance-Vought  Corp. 

Post  Office  Box  5907 
Dallas,  Texas 

Library  1-6310/3L-2884 

National  Science  Foundation 
1951  Constitution  Avenue,  N.  W, 

Washington  25,  D.  C. 

Attn:  Engineering  Sciences  Division 

New  York  University 
University  Heights 
New  York  53,  New  York 

Attn:  Department  of  Aeronautical  Engineering 
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New  York  University 
25  Waverly  Place 
New  York,  New  York 

Attn:  Library,  Institute  of  Math.  Sciences 

NORAIR 

A  Division  of  Northrop  Corporation 
Hawthorne,  California 
Attn:  Library 

Northrop  Aircraft, ,  Inc. 

Hawthorne,  California 
Attn:  Library 

Gas  Dynamics  Laboratory 
Technological  Institute 
Northwestern  University 
Evanston,  Illinois 
Attn :  Library 

Pennsylvania  State  University 
University  Park,  Pennsylvania 

Attn:  Library,  Dept,  of  Aero.  Engineering 

The  Ramo-Wooldridge  Corporation 

8820  Bellanca  Avenue 

Los  Angeles  45,  California 

Gifts  and  Exchanges 
Fondren  Library 
Rice  Institute 
P.  0.  Box  1892 
Houston,  Texas 

University  of  Southern  California 
Engineering  Center 
Los  Angeles  7,  California 
Attn:  Librarian 

The  Editor 

Battelle  Technical  Review 
Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  Ohio 

Douglas  Aircraft  Company,  Inc. 

Long  Beach,  California 
Attn:  Library 
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FluiDyne  Engineering  Corporation 
5740  Wayzata  Boulevard 
Golden  Valley 
Minneapolis ,  Minnesota 

Grumman  Aircraft  Engineering  Corporation 
Bethpage,  Long  Island,  New  York 

Lockheed  Missiles  and  Space  Company 
P.  0,  Box  551 
Burbank,  California 
Attn:  Library 

Marquardt  Aircraft  Corporation 

7801  Havenhurst 

Van  Nuys,  California 

Martin  Company 
Denver,  Colorado 

Martin  Company 
Orlando,  Florida 
Attn :  J.  Mayer 

Mississippi  State  College 

Engineering  and  Industrial  Research  Station 

Aerophysics  Department 

P.  0.  Eix  248 

State  College,  Mississippi 


Lock.ueed  Missiles  and  Space  Company 
3251  Hanover  Street 
Palo  Alto,  California 
Attn:  Library 

General  Electric  Company 
Research  Laboratory 
Schenectady,  New  York 

Attn:  Dr.  H,  T.  Nagamatsu 
Attn :  Library 

Fluid  Dynamics  Laboratory 
Mechanical  Engineering  Department 
Stevens  Institute  of  Technology 
Hoboken,  New  Jersey 

Attn:  Dr.  R.  H.  Page,  Director 
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Department  of  Mechanical  Engineering 
University  of  Arizona 
Tucson,  Arizona 

Attn:  Dr.  E.  K.  Parks 

Vitro  Laboratories 
200  Pleasant  Valley  Way 
West  Orange,  New  Jersey 

Department  of  Aeronautical  Engineering 
University  of  Washington 
Seattle,  Washington 

Attn:  Prof.  R.  E.  Street 
Attn:  Library 

American  Institute  of  Aeronautics  &  Astronautics 
1290  Avenue  of  the  Americas 
New  York,  New  York 

Attn:  Managing  Editor 
Attn:  Library 

Department  of  Aeronautics 
United  States  Air  Force  Academy 
Colorado 

MHD  Research,  Inc. 

Newport  Beach,  California 
Attn:  Technical  Director 

University  of  Alabama 
College  of,  Engineering 
University,  Alabama 

Attn:  Head,  Dept,  of  Aeronautical 
Engineering 

ARDE  Associates 
100  W.  Century  Road 
Pararaus ,  New  Jersey 

Attn:  Mr.  Edward  Cooperraan 

Aeronautical  Research  Associates  of  Princeton 
50  Washington  Road 
Princeton,  New  Jersey 

Attn:  Dr,  C.  duP.  Donaldson,  President 
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Daniel  Guggenheim  School  of  Aeronautics 
Georgia  Institute  of  Technology 
Atlanta,  Georgia 

Attn:  Prof.  A.  L.  Ducoffe 

University  of  Cincinnati 
Cincinnati,  Ohio 

Attn:  Prof.  R.  P.  Harrington,  Head 

Dept,  of  Aeronautical  Engineering 
Attn:  Prof.  Ting  Yi  Li,  Aerospace  Engineering  Dept. 

Virginia  Polytechnic  Institute 
Dept,  of  Aerospace  Engineering 
Blacksburg,  Virginia 
Attn:  Dr.  R.  T.  Keefe 
Attn:  Dr.  J.  B.  Eades,  Jr. 

Attn:  Library 

IBM  Federal  System  Division 
7220  Wisconsin  Avenue 
Bethesda,  Maryland 

Attn:  Dr.  I.  Korobkin 

Superintendent 

U.  S.  Naval  Postgraduate  School 
Monterey,  California 

Attn:  Technical  Reports  Section  Library 

National  Bureau  of  Standards 
Washington  25,  D.  C. 

Attn:  Chief,  Fluid  Mechanics  Section 

North  Carolina  State  College 
Raleigh,  North  Carolina 

Attn:  Division  of  Engineering  Research 
Technical  Library 

Defense  Research  Corporation 
P.  0.  Box  No.  3587 
Santa  Barbara,  California 
Attn:  Dr.  J.  A.  Laurmann 

Aerojet -General  Corporation 
6352  North  Irvindale  Avenue 
Box  296 

Azusa,  California 
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Apollo  -  DDCS 
General  Electric  Company 
A&E  Building,  Room  204 
Daytona  Beach ,  Florida 
Attn:  Dave  Hovis 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis ,  Minnesota 

Attn:  Prof.  J.  D.  Akerman 

Guggenheim  Laboratory 
Stanford  University 
Stanford ,  California 

Attn:  Prof.  D.  Bershader,  Department  of  Aero. 
Engineering 

Space  Technology  Laboratory,  Inc. 

1  Space  Park 

Redondo  Beach,  California  90200 

Attn:  STL  Tech.  Lib.  Doc.  Acquisitions 

University  of  Illinois 

Department  of  Aeronautical  and  Astronautical  Engineering 
Urbana,  Illinois 

Attn:  Prof.  H.  S.  Stilwell 

Armour  Research  Foundation 
Illinois  Institute  of  Technology 
10  West  35th  Street 
Chicago,  Illinois 

Attn:  Dr.  L.  N.  Wilson 

Institute  of  the  Aeronautical  Sciences 
Pacific  Aeronautical  Library 
7600  Beverly  Boulevard 
Los  Angeles ,  California 

University  of  California 
Department  of  Mathematics 
Los  Angeles ,  California 
Attn:  Prof.  A.  Robinson 

Louisiana  State  University 
Department  of  Aeronautical  Engineering 
College  of  Engineering 
Baton  Rouge,  Louisiana 
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Mathematical  Reviews 
American  Mathematical  Society 
80  Waterman  Street 
Providence,  Rhode  Island 

Stanford  University 

Department  of  Aeronautical  Engineering 
Stanford,  California 
Attn :  Library 

University  of  California 
Aeronautical  Sciences  Laboratory 
Richmond  Field  Station 
1301  South  46th  Street 
Richmond ,  California 

University  of  Denver 

Department  of  Aeronautical  Engineering 

Denver  10,  Colorado 

University  of  Chicago 
Laboratories  for  Applied  Sciences 
Museum  of  Science  and  Industry 
Chicago,  Illinois 
Attn:  Librarian 

University  of  Colorado 

Department  of  Aeronautical  Engineering 

Boulder,  Colorado 

University  of  Illinois 
Aeronautical  Department 
Chamnaign,  Illinois 

University  of  Kentucky 

Department  of  Aeronautical  Engineering 

College  of  Engineering 

Lexington,  Kentucky 

University  of  Toledo 

Department  of  Aeronautical  Engineering 
Research  Foundation 
Toledo,  Ohio 
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Aerospace  Corporation 

P.  0.  Box  95085 

Los  Angeles ,  California 

Attn:  Advanced  Propulsion  &  Fluid  Mechanics  Department 
Attn:  Gas  Dynamics  Department 

Boeing  Scientific  Research  Laboratory 
P.  0.  Box  3981 
Seattle,  Washington 

Attn:  Dr.  A.  K.  Sreekanth 
Attn:  G.  J,  Appenheimer 

Vidya,  Inc. 

2626  Hanover 
Palo  Alto,  California 
Attn:  Mr.  J.  R.  Stalder 
Attn :  Library 

General  Electric  Company 

FPD  Technical  Information  Center  F-22 

Cincinnati,  Ohio 

Northwestern  University 
Technological  Institute 
Evanston,  Illinois 

Attn:  Department  of  Mechanical  Engineering 

Harvard  University 
Cambridge,  Massachusetts 

Attn:  Prof,  of  Engineering  Sciences  &  Applied  Physics 
Attn :  Library 

University  of  Wisconsin 
P.  0.  Box  2127 
Madison,  Wisconsin 

Attn:  Prof.  J.  0.  Hirschf elder 

Dr.  Antonio  Ferri,  Director 
Guggenheim  Aerospace  Laboratories 
New  York  University 
181st  St.  and  University  Ave. 

Bronx,  New  York 
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Department  of  Aerospace  &  Mechanical 
Engineering  Sciences 
University  of  California-San  Diego 
La  Jolla,  California  92037 
Attn:  Dr.  P.  A.  Libby 

University  of  New  Mexico 
Albuquerque,  New  Mexico 

Attn:  Dr.  Theodore  Sparks 

Notre  Dame  University 
Southbend,  Indiana 

Attn:  Dr.  John  D.  Nicolaides 

Department  of  Aerospace  Engineering 

Dr.  John  Laufer,  Chairman 

University  of  Southern  California 

Graduate  Department  of  Aerospace  Studies 

University  Park 

Los  Angeles,  California  90007 
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